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ABSTRACT 

We present a compositional analysis of the 10 fim silicate spectra for brown 
dwarf disks in the Taurus and Upper Scorpius (UppSco) star-forming regions, 
using archival Spitzer/IRS observations. A variety in the silicate features is ob- 
served, ranging from a narrow profile with a peak at 9.8 fim, to nearly flat, 
low-contrast features. For most objects, we find nearly equal fractions for the 
large-grain and crystalline mass fractions, indicating both processes to be active 
in these disks. The median crystalline mass fraction for the Taurus brown dwarfs 
is found to be 20%, a factor of ~2 higher than the median reported for the higher 
mass stars in Taurus. The large-grain mass fractions are found to increase with 
an increasing strength in the X-ray emission, while the opposite trend is observed 
for the crystalline mass fractions. A small 5% of the Taurus brown dwarfs are still 
found to be dominated by pristine ISM-like dust, with an amorphous sub-micron 
grain mass fraction of ~87%. For 15% of the objects, we find a negligible large- 
grain mass fraction, but a >60% small amorphous silicate fraction. These may 
be the cases where substantial grain growth and dust sedimentation has occurred 
in the disks, resulting in a high fraction of amorphous sub-micron grains in the 
disk surface. Among the UppSco brown dwarfs, only usdl61939 has a S/N high 
enough to properly model its silicate spectrum. We find a 74% small amorphous 
grain and a ~26% crystalline mass fraction for this object. 

Subject headings: stars: abundances — circumstellar matter — stars: low-mass, 
brown dwarfs 



1. Introduction 

Several surveys in recent years in the near- and mid-infrared have confirmed the presence 
of circum(sub)stellar disks around brown dwarfs (e.g., Muench et al. 2001; Liu et al. 2003; 
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Luhman et al. 2005), and indicate a range in disk properties, similar to T Tauri disks. 
The disk fractions for stellar and sub-stellar members in young star-forming regions are 
found to be comparable (e.g., Liu et al. 2003; Luhman et al. 2005), and may even exceed 
the fractions found for higher mass T Tauri stars (e.g., Scholz et al. 2007; Riaz & Gizis 
2008). A key question in the study of disk evolution is the extent of dust processing that 
has occurred in the disk, as indicated by processes such as grain growth, crystallization 
and dust settling. In this regard, a detailed compositional analysis of the 10 /im silicate 
emission feature has been immensely valuable in gaining insight into such processes, that 
may eventually lead to planet formation. The most dominant dust species in the circumstellar 
material are oxygen- rich silicates, that are either olivines (Mg2j;Fe2(i_a;)Si04), ranging from 
fayalite {x = 0) to forstcritc (,t = 1), or pyroxenes (Mg2:Fe(i_2.)Si03), ranging from ferrosillite 
{x = 0) to cnstatitc {x = 1) (e.g., van Bockel ct al. 2005; Kcsslcr-Silacci et al. 2005). 
Silicate dust in the interstellar medium (ISM) is dominated by amorphous olivine, the 10 
/xm spectrum of which shows a narrow peak at 9.8 fim (e.g., Bouwman et al. 2001). In 
more evolved systems such as the comet Hale-Bopp or the debris disk around /3 Pictoris, 
the silicate spectrum is dominated by crystalline enstatite and forsterite (e.g., Crovisier et 
al. 1997; Bouwman et al. 2001; Okamoto et al. 2004) that show prominent emissivity 
peaks near 9.3 and 11.3 /im, respectively. An evolutionary sequence from ISM-like dust to 
cometary dust thus involves a composition change from amorphous to crystalline silicates, 
indicated by a shift in the peak wavelength of the 10/im spectrum from ~9.8 to 11.3//m 
(e.g., Bouwman et al. 2001; Kcssler-Silacci et al. 2005). This is also accompanied by a 
broadening and flattening of the observed spectrum, as smaller grains coagulate into larger 
particles and settle to the optically thick disk midplane. Detailed studies on the extent 
of grain processing in Herbig Ae/Be systems have indicated a lack of ISM-like dust, and 
a correlation between the grain growth and crystallization processes (e.g., Bouwman et al. 
2001; van Boekel et al. 2005). Bouwman et al. (2001) also report on an increase in the 
abundance of silica (Si02) with increasing crystallinity, another important species formed 
when amorphous silicates are thermally annealed to form forsterite (e.g., Fabian et al. 2000). 
Among lower mass (M* < l^o) T Tauri systems, a trend of decreasing small and crystalline 
grain abundances is observed with age, and both the large-grain and crystalline mass fractions 
arc found to increase with the degree of dust settling in the disk (e.g., Sargent et al. 2009; 
Watson ct al. 2009; Bouwman ct al. 2008). Honda et al. (2006) report crystallization of 
5%-20% amorphous sihcate grains in T Tauri systems, regardless of the stellar age. Dust 
compositional analyses of sub-stellar objects indicate a wide range in the large-grain and 
crystalline mass fractions, with a degree of crystallinity as high as that observed among 
Herbig Ae/Be stars (e.g., Apai et al. 2005). A more flattened disk structure is observed 
for the disks around cool stars, suggesting a higher degree of dust processing in these disks 
compared to higher mass stars (e.g., Apai et al. 2005; Kessler-Silacci et al. 2006; Pascucci 
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et al. 2009). 

We present Spitzer Infrsned Spectrograph (IRS) archival observations of brown dwarfs in 
the Taurus-Auriga (~l-2 Myr; Kenyon & Hartmann 1995) and the Upper Scorpius (UppSco) 
(~5 Myr, Preibisch & Zinnecker 1999) star-forming regions. Combining our work with the 
results from Apai et al. (2005) in the Chamaeleon I (Cha I) star-forming region, we compare 
the mass fractions for crystalline and large silicates among brown dwarfs at different ages, 
as well as with higher mass stars, and discuss any possible correlations observed. 

2. Targets 

Our criteria for selecting targets in Taurus was to search for spectroscopically confirmed 
members with spectral types (SpT) later than M6, making these more likely to be sub- 
stellar objects (e.g., Luhman 2004), and for which excess emission has been detected in 
the Spitzer/IRAC bands, thus confirming the presence of disks around them. We therefore 
mainly searched for objects discussed in the works by Luhman et al. (2006) and Guieu et 
al. (2007). We were able to obtain publicly released Spitzer/IRS observations for 20 objects. 
Five of these do not show any emission in the lO/xm silicate feature and have flat spectra, 
i.e. the emission levels are consistent with the underlying continuum. The 15 objects with 
emission features are listed in Table 1. Fig. [T] shows the normalized continuum-subtracted 
spectra in units of [{F^ — Fc)/Fc] for these objects, along with the 5 flat spectra shown in 
the bottom panel. The method for normalizing the silicate feature to the continuum level is 
discussed in §4.21 In UppSco, Scholz et al. (2007) have reported the presence of a solid-state 
10 /xm emission feature in only 3 out of the 13 brown dwarfs with confirmed disk excess. 
These are also listed in Table 1, with the continuum-subtracted spectra shown in Fig. [T] 
The Spitzer PIDs for these observations are 30540, 2, 248, and 20435. We note that four of 
these objects have SpT of M5.5 (with an uncertainty of ±0.25), but since these have been 
discussed earlier in the literature (e.g., Luhman et al. 2006; Luhman 2006; Scholz et al. 
2007) as brown dwarf members, we have included them in our study. 

The IRS spectra presented here were obtained in the low spectral resolution {X/6X ~ 
90) Short-Low [SL] module that covers 5.2-14.5 /im in 2 orders. The observations were 
obtained at two nod positions along each slit. The Basic Calibrated Data (BCD) were 
produced by the S15 pipeline at the Spitzer Science Center (SSC). The sky background 
was removed from each spectrum by subtracting observations taken in the same module, 
but with the target in the other order. The spectra were extracted and calibrated using 
the Spitzer IRS Custom Extraction (SPICE) software provided by the SSC. The extraction 
was performed using a variable-width column extraction that scales with the width of the 
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wavelength-dependent point spread function. The spectra at the two nod positions were 
then averaged to obtain a single spectrum. The signal-to-noise (S/N) ratio for the spectra 
varies strongly with the wavelength, and the noise levels are found to be higher at the red 
end. For the spectra presented here, the S/N varies between ~8 and >50; the more noisy 
spectra such as 2M04554801 and uscoll2 have S/N values that are quite low (<10), while 
it is between ~30 and 40 for MHO 5 and 2M04141188. 



3. Method 

Five dust species were considered in order to obtain a compositional fit to the observed 
spectra: amorphous olivine and pyroxene, crystalline enstatite, crystalline forsterite, and 
silica. The characteristics of these are outlined in Table 2, and the spectral profiles plotted in 
Fig. [21 The bulk densities listed in Table 2 have been obtained from Bouwman et al. (2001). 
The amorphous grains show little dependence on shape (e.g., Bouwman et al. 2001), and 
thus can be assumed to be homogeneous and spherical so that the standard Mie theory can 
be applied to determine their spectroscopic properties. The properties of crystalline silicates 
however are sensitive to the grain shape which makes it important to consider shapes other 
than homogeneous spheres. Min et al. (2003) have studied the absorption and scattering 
properties of particles that are inhomogeneous in structure and composition, based on the 
assumption that the characteristics of such irregularly shaped particles can be simulated by 
the average properties for a distribution of shapes, such as ellipsoids, spheroids or hollow 
spheres. We considered both the distribution of hollow spheres (DHS) and the continuous 
distribution of ellipsoids (CDE) from Min et al. (2003), and found a better match both in 
terms of the location and width of the spectral profiles using the CDE routine. Therefore the 
absorption efficiencies of crystalline silicates (enstatite and forsterite) and silica have been 
calculated using ellipsoidal grains. Since the CDE method has been used, the grains are 
treated in the "Rayleigh limit" , i.e. the grain sizes are assumed to be much smaller than the 
wavelength of radiation. Our model is therefore biased towards sub- micron (0.1 fim) sized 
crystalline grains, and the effects of grain growth could not be studied for these species. 
For the case of amorphous silicates, two grain sizes of 0.1 and 2.0 /xm have been used. As 
discussed in Bouwman et al. (2001), the optical properties of dust grains with a range in 
sizes between 0.01 and 5 fim can be characterized by two typical grain sizes: 0.1 /zm for 
"small" grains with sizes < 1 fim, and 2.0 /im for "large" grains with sizes >1 fim. 

To construct a model spectrum, we have followed the method outlined in Honda et al. 
(2003), wherein a power-law source function is assumed and the model flux Fx is the sum 
of the emission from a featureless power-law continuum and the emission from the different 
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dust species. The model spectrum can be written as 
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+ (OlQo.lao + 0'2Q2.Qao + ^sQo.lap + 0'iQ2.Qap 




where Qo.iao and Q2.oao are the absorption efficiencies for 0.1 and 2.0 /im amorphous olivine, 
Qo.iap and Q2.oap are the absorption efficiencies for 0.1 and 2.0 /im amorphous pyroxene, 
and Qens, Qfors, and Qsu are the absorption efficiencies for enstatite, forsterite and sihca, 
respectively. There are 10 free parameters; eight multiplicative factors Oo-a/ in 10"^"^ W 
m~^, and two spectral indices of the source function, m and n. To fit the spectrum, we used 
the x^-minimization method outlined in van Boekel et al. (2005). The best-fit values thus 
obtained for the multiplicative factors are listed in Table 1, along with the derived mass 
fractions. Also listed is the reduced-x^ value (x^/dof; dof=651) of the fit. The uncertainties 
for the fit parameters were obtained by considering the errors on the observed fiuxes, as well 
as the S/N ratio. We used the Monte Carlo method for error estimation as outlined in van 
Boekel et al. (2005), and discussed further in Juhasz et al. (2009). In this method, random 
Gaussian noise is added to the spectrum, with 1000 synthetic spectra generated and fitted 
at each noise level. We considered four different noise levels with a p^^rror ^ -pobserved q 
0.1, 0.01 and 0.001, corresponding to a S/N ratio of 5, 10, 100 and 1000, respectively. The 
mean of the resulting distribution of all fit parameters then corresponds to the best-fit value, 
and the errors are derived from the standard deviation. 

In order to validate our method, we modeled some known low-mass young objects and 
compared our results with those published in the literature. Fig. [3] shows the observed and 
modeled spectra for two such objects at SpT~M6, V410 Anon 13 in Taurus and SSTc2d 
J161159.9-382337 (SST-Lup3-1) in the Lupus III dark cloud. V410 Anon 13 shows a high 
abundancy of crystalline silicates but a negligible large grain mass fraction (Sargent et al. 
2009). We find a small amorphous, large amorphous and crystalline silicate mass fraction 
of 40.2%!^, 0.002% ± 0.1%, and 44.3%!^ for this object, which is consistent within the 
uncertainties with the fractions of 46.3%ibl7.2%, 10.4%±20.1%, and 37.2%±11.3% found 
by Sargent et al. (2009). Sargent et al. (2006) had earlier reported a 0% large-grain 
mass fraction for V410 Anon 13. For SST-Lup3-1, Merm et al. (2007) have reported a 
~33% crystalline mass fraction, and an equal mix of small and large grains from the hot 
component (T ~ 300 K; high-continuum case) of their compositional fit. We find a small 
and large amorphous grain mass fraction of 40.3%l4 and 37.4:%~^l for this object, with a 
crystalline (sub-micron) mass fraction of 28.3%l3, similar to the one reported by Merin et 
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al. (2007). Our derived crystalline mass fractions for some of the higher mass T Tauri stars 
from Sargent et al. (2009) work, as well as the Cha I brown dwarfs from Apai et al. (2005), 
were also found to be consistent within the uncertainties (better than 8%) with the fractions 
reported by these authors. We note that since our model does not include opacities of large 
crystalline grains, it underestimates the contribution to the total crystalline mass fraction 
due to large-sized grains. This makes a direct comparison of our 28.3% crystalline mass 
fraction for SST-Lup3-1 with the 33% reported by Merm et al. difficult, since they find the 
crystalline silicates to be a mixture of small and large grains (in the hot component). We 
refer to the discussion in Sargent et al. (2009) on the opacities of large crystalline grains. 
These authors have presented a compositional analysis of T Tauri stars in Taurus, and have 
also considered only sub-micron crystalline silicates in their model. As explained by these 
authors, large grains should be considered as heterogenous aggregates, made up of sub-micron 
sized amorphous silicates, forsterite or enstatite, with the highest abundance of amorphous 
grains. The opacities of the forsterite components of such large heterogenous grains have been 
shown by Min et al. (2008) to resemble the opacity of small homogenous forsterite grains, 
whereas the opacity of the amorphous component of the aggregate is found to resemble that 
of large amorphous homogenous grains, with a size comparable to that of the heterogenous 
aggregate. Therefore modeling the crystalline component of a large heterogenous aggregate 
with small crystalline grains should not result in largely erroneous mass fractions for the 
crystalline silicates. This, as Sargent et al. (2009) indicate, is consistent with the findings 
of Bouwman et al. (2008), who have considered three different grain sizes of 0.1, 1.5 and 6.0 
Hm in their models, but find the typical sizes from the model fits to be sub-micron for the 
crystalline silicates, and 6 fim for the amorphous grains. 



4. Results 

Fig. Hlshows the models fits obtained for the Taurus and UppSco objects, along with the 
contribution from the individual components. Among the Taurus brown dwarfs, 2M04141760 
shows a narrower feature compared to other objects, with a peak in emission near 9.8 fim 
(Fig. [1]). It exhibits the least signs of dust processing, with a small amorphous silicate mass 
fraction of ~87±10% and large grain and crystalline mass fractions of <10%. This object 
is similar to the case of Cha Ha 1 that shows a smooth narrow feature with a peak at 9.8 
/im (Apai et al. 2005). CFHT-BD-Tau 8, V410 X-ray 6, CFHT-BD-Tau 6, 2M04400067 and 
2M04230607 also show peak positions near 9.8 or 9.3 /im, and have small amorphous silicate 
mass fractions of >;50%. V410 X-ray 6 also has the lowest crystalline mass fraction (0.002%) 
in the sample. The object 2M04141188 shows a clear peak in emission near 11.3 fim (Fig. [1]) 
and a broad, fiat-peaked feature; we find a ~20% crystalline mass fraction for this object 
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and roughly equal fractions between 30 and 40% for the small and large amorphous grains. 
MHO 5, 2M04414825 and CFHT-BD-Tau 9 also exhibit broad, fiat features. MHO 5 and 
2M04414825 show the strongest signs of grain growth in their disks, with large grain mass 
fractions of 64% and 73%, respectively, while a higher degree of crystallinity (32% by mass) 
must be responsible for the flattening observed in CFHT-BD-Tau 9. High crystalline mass 
fractions between 40% and 50% are found for 2M04290068 and CFHT-BD-Tau 20. The case 
for CFHT-BD-Tau 20 however is dubious given its noisy spectrum. It shows the weakest 
feature among the Taurus members, and has a negligible large grain mass fraction. Among 
the UppSco objects, only usdl61939 has a S/N high enough to obtain good estimates on 
the mass fractions. Modeling of its silicate spectrum indicates that the disk may still be 
dominated by small amorphous grains (74% by mass). A high crystalline mass fraction is 
also found for uscoll2 (~82%), but similar to CFHT-BD-Tau 20, the spectrum is too weak 
and noisy which makes the derived mass fractions highly uncertain. 



4.1. Additional Model Fits 



For the three UppSco objects and the Taurus brown dwarf CFHT-BD-Tau 20, we con- 
structed additional model spectra considering only amorphous silicates (both small- and 
large-sized) in the model, and excluding all crystalline grains (Fig. Ek). This case '(a)' there- 
fore highlights the importance of crystalline silicates in these disks, particularly in the case 
of USC0II2 and CFHT-BD-Tau 20 that are found to have high crystalline mass fractions 
(^50%), despite the absence of any clear crystalline features in their observed spectra. In 
addition, these objects show a negligible amount of large grains in their disks (mass frac- 
tion <0.1%). We therefore constructed another set '(b)' of model spectra considering only 
large-sized amorphous grains (Fig. third panel). This set (b) was also constructed for 
the objects CFHT-BD-Tau 8, 2M04230607 and 2M04400067 that have large grain mass frac- 
tions between 0.002% and 0.003% (Fig. [Sjo). We have listed for comparison in Table 1 the 
reduced-x^ values for these additional model-fits. The first panels in these figures show the 
best-fits obtained (based on the lowest reduced-x^ value) using all five species. The fits for 
cases (a) and (b) for uscoll2 result in an increase in the reduced-x^ value of only 0.03 and 
0.08, respectively, indicating that all three could fit the observed spectrum well. The crys- 
talline fraction for this object therefore could be as low as 0%, and the large grain fraction 
as high as 100%. For the case of usdl60958, the reduced-x^ value is the same for the best-fit 
and the case (a) fit, indicating crystalline silicates to be insignificant in this disk. The fit for 
case (b) though has a larger reduced-x^ value of 1.7, and shows that some fraction of small 
grains is required to fit the observed spectrum. For the case of usdl61939, the difference in 
the reduced-x^ values for the best-fit and case (a) fit is again very small (0.04), though the 
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case (b) model is not a good fit. Inspecting the fits indicates that some fraction of crystalhne 
sihcates should be included to better fit the spectrum near 11.3 fim. The small amorphous 
fraction for this object could therefore lie between ~74 and 100%, and the crystalline frac- 
tion between and 26% (Table 1). CFHT-BD-Tau 20 has the flattest feature among Taurus 
objects. Fits for the cases (a) and (b) result in an increase in the reduced-x^ value of 0.09 and 
0.23, respectively. Inspecting the fits closely again indicates that some crystalhne fraction 
is required to better fit the observed spectrum near 11.3 /im, though the feature is so fiat 
that large grains alone could fit the spectrum. In the case of CFHT-BD-Tau 8, 2M04230607 
and 2M04400067 (Fig. [5]d), the case (b) fits result in an excess emission around ~12 fim, 
and are also bad fits near ~9 fim for 2M04230607 and 2M04400067. The differences in the 
reduced-x^ values are between ~0.3 and 0.7. A model with large grains only therefore is 
inadequate to fit these spectra, and the negligible large grain mass fractions as deduced from 
the best-fits must then be good estimates for these objects. 

4.2. Shape vs. Strength of the Silicate Features: Dependence on Grain 
Grovv^th, Crystallization and Dust Settling 

In the discussion that follows, we have excluded the objects uscoll2, usdl60958 and 
CFHT-BD-Tau 20 due to the large uncertainties in their model-fits and the derived mass 
fractions. Figure [H^ compares the 'shape' and 'strength' in the 10 fim silicate feature for 
Taurus and UppSco brown dwarfs with sun-like T Tauri stars in Taurus (Pascucci et al. 
2009) and Herbig Ae/Be stars (van Boekel et al. 2005). The feature strength can be 
estimated from the peak normalized fiux, Fp^ak, above the continuum, while the shape can 
be estimated by the parameter Fn ^/Fg ^, the ratio of the normalized fiuxes at 11.3 and 9.8 
fim (e.g., Bouwman et al. 2001). In order to normalize the silicate feature to the continuum, 
we first approximated the underlying continuum by connecting the two end points of the 
spectrum, i.e. the observed fiuxes at 8 and 13 fim (averaging over 0.2 fim at the end points; 
Kessler-Silacci et al. 2005). The normalized silicate spectrum is then defined in units of 
[{Fi, — Fc)/Fc], where F^, and F^ are the observed and the continuum normalized fiuxes, 
respectively. The continuum-subtracted normalized spectra are shown in Fig. [H and the 
values for Fpeak and -Fn.s/Fg.g plotted in Fig. [6^ have been calculated from these normalized 
spectra. Typical uncertainties in Fpeak are ±0.1 for the brown dwarfs, ±0.3 for the T Tauri 
stars, and ±0.6 for the Herbig Ae/Be stars, while the uncertainties in the Fn ^/Fg ^ values 
are ±0.1 for the brown dwarfs and T Tauri stars, and ±0.15 for the Herbig Ae/Be stars, van 
Boekel et al. (2003) first noted a correlation between the shape and strength of the emission 
features for a sample of Herbig Ae/Be stars. The correlation was such that spectra with a 
larger value for the 11.3 to 9.8 fiux ratio were found to have a lower peak-over-continuum fiux. 
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This was explained by the extent of dust processing in the disk; spectra that are dominated 
by more processed dust due to grain growth and/or crystaUization processes were found to 
be flatter with a peak wavelength close to 11.3 /xm, while spectra dominated by unprocessed 
ISM-like dust showed more peaked, narrower profiles with a peak wavelength closer to 9.8 
/im. Similar correlation have been observed among T Tauri stars (e.g., Pryzgodda et al. 
2003; Kessler-Silacci et al. 2005; Pascucci et al. 2009) and sub-stellar objects (Apai et al. 
2005; Pascucci et al. 2009). In a strength vs. shape diagram then, the narrow, peaked 
profiles would occupy the space near the bottom right corner, while with increasing dust 
processing the feature would move towards the upper left region (e.g., Kessler-Silacci et al. 
2006). Apai et al. (2005) further noted a 'reversal' in this trend for Cha I brown dwarfs, such 
that for highly processed dust the feature strength may continue to increase with increasing 
crystallinity. 

Fig. [6^1 shows the Taurus and UppSco brown dwarfs clustered around Fp^ak ~ 1.3 and 
F11.3/F9.8 ~ 1.0, with the features being weaker in these objects compared to most other 
T Tauri and Herbig Ae/Be stars. The brown dwarf 2M04141760 lies to the right of all 
other objects in the sample, as expected given the stronger, more peaked feature observed 
for this object (Fig. [T]). We do not observe any of the extreme cases as in Cha I with 
strong features for highly processed grains. The presence of weaker features in our objects 
is consistent with the findings from previous works (e.g., Apai et al. 2005; Kessler-Silacci et 
al. 2006; Sargent et al. 2009; Pascucci et al. 2009): disks around cool stars exhibit fiatter 
silicate emission features, indicating a higher degree of dust processing compared to the disks 
around higher mass stars. Kessler-Silacci et al. (2007) have explained this by the differences 
in the location of the silicate emission zone for stars of different luminosities. The 10 fim 
feature probes smaller radii (-R10 < 0.001-0.1 AU) in disks around brown dwarfs than in 
disks around T Tauri stars {Rio > 0.1-3 AU). In a recent compositional analysis of T Tauri 
disks in Taurus, Sargent et al. (2009) have reported higher abundances of large grains in the 
"inner" disk regions probing radii of ~0.6 AU, than in the "outer" regions that probe larger 
disk radii around 10 AU, suggesting that the grain growth process occurs more rapidly at 
smaller radii. For the crystalline silicates, the abundances were found to be similar in both 
regions. Thus due to the smaller radii probed by the 10 /xm feature around brown dwarfs, 
stronger signatures of dust processing would be observed, that may or may not represent the 
properties of the bulk silicate dust in the disk (Kessler-Silacci et al. 2007). 

In Figs. [6]b-e, we investigate whether the observed fiattening and the shift in the peak 
position of the silicate features is due to an increase in the grain sizes or a higher degree of 
crystallinity in the disk. No clear rise in the grain growth or the crystallinity levels is observed 
towards decreasing Fp^ak values (Figs. [6]b-c), and the fractions are largely dispersed for Fpeak 
between ~1.5 and 1.3. Thus neither processes show any correlation with the strength in the 
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silicate feature. One exception is 2M04141760 that shows the narrowest sihcate profile, 
consistent with a small ~8% and 5% large and crystalline mass fraction, respectively. In 
Figs. E] d-e, both the large grain and crystalline mass fractions are found to increase as 
the 11.3 fim fiux increases relative to the 9.8 /xm fiux. Thus the feature shape seems to 
be correlated with both processes. This is contrary to the results seen among higher mass 
objects in Taurus, as Watson et al. (2009) report both the shape and strength to be strongly 
linked to the crystallinity in these disks, but no strong correlation is found for the large grain 
abundances. Watson et al. have discussed these spectral indices not to be good tracers of the 
large grain mass fractions, since the shape-strength correlation is mainly due to the trend 
of sedimentation and crystallinity, and is thus a good tracer of grain growth only in objects 
with small crystalline mass fractions. The correlations observed in Figs. [6] d-e nevertheless 
suggest both the grain growth and crystallization processes to be active in these brown dwarf 
disks, as both seem to affect the shape of the observed features. 

Fig. shows the extent of flaring in the brown dwarf disks in our sample. Here, we 
plot the flux densities at 8 and 13 /im, normalized to the flux at 8 /xm. These wavelengths 
mark the end points of the 10 /zm silicate emission feature, and have been used to determine 
the continuum in the feature. As discussed in Apai et al. (2005), the continuum of the 
silicate spectrum is deflned by the disk geometry, and therefore the normalized flux density 
at 13 fim will be higher for a highly flared disk, and will decrease as the disk geometry 
gets flatter. Disks dominated by pristine ISM-like dust show flared structures, while flatter 
structures are observed for the ones that have gone through substantial grain processing. 
The flattening is mainly caused by gravitational settling of dust grains from the optically 
thin upper atmosphere to the optically thick disk midplane. The processing of dust into 
crystalline silicates as well as growth to larger sizes both effect the vertical structure of the 
disk. Sargent et al. (2009) and Watson et al. (2009) have reported an increase in the 
mass fractions for both crystalline and large silicates as the degree to which the disks are 
sedimented increases. Among the Taurus brown dwarfs, 2M04141760 shows the most flared 
geometry, consistent with an amorphous sub-micron grain mass fraction of ~87%. The extent 
of flaring in this disk is similar to Cha Ha 1 that is also dominated by unprocessed dust grains 
(<10% crystalline mass fraction; Apai et al. 2005). Most other Taurus brown dwarfs, along 
with the UppSco object usdl61939, show intermediate flaring, with the flattest structures 
observed for MHO 5 and GM Tau. Comparing Figs. [7^ and b, MHO 5 flat geometry is 
consistent with a high large-grain mass fraction of ~64%. However, 2M04414825 that shows 
the strongest signs of grain growth in the sample has a more flared disk than MHO 5 or GM 
Tau. Similar is the case for 2M04290068 that has the highest crystalline mass fraction in the 
sample but shows intermediate flaring. The brown dwarf GM Tau has the flattest structure 
but shows nearly equal fractions for large and crystalline silicates. No strong correlation is 
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therefore observed for these fractions with the extent of dust setthng in the disk. 

The sedimentation process however is effected by the level of turbulence in the disk. 
In disks with higher accretion rates, turbulent mixing may occur alongside grain growth 
and thus prolong the timescales over which large grains settle to the optically thick disk 
midplane (e.g., DuUemond & Dominik 2004). As discussed in DuUemond & Dominik (2008), 
with decreasing level of turbulence, the level of sedimentation would increase, resulting in an 
under-abundance of large grains in the upper disk layers. In Fig. [Tt, we observe the opposite 
for the case of MHO 5 and 2M04414825 that have weaker accretion rates than the rest but 
show high large-grain mass fractions of ~60-70%. There are however two more objects in 
Fig. [Tt that show small fractions of ~20% at similar accretion rates. The more actively 
accreting disk 2M04141760 shows a flared geometry, but has a small 8% large-grain mass 
fraction. In Fig. [TJi, no clear correlation is observed between crystallization and the disk 
mass accretion rate, as similar levels of crystallinity are found for logM between ~-9 and -11 
Mq yr^^. Among higher mass T Tauri stars in Taurus, while Sargent et al. (2009) have noted 
higher abundances for the crystalline and large silicates in the more settled disks, Pascucci 
et al. (2009) have found similar large-grain fractions at all accretion rates, thus suggesting 
that turbulence in the disk may have little or no effect on the dust processing mechanisms. 
A larger number of data points would be valuable to investigate any dependencies on the 
mass accretion rate. 

Returning to Fig. [Tb, most objects show roughly equal fractions for the large and crys- 
talline grains, again indicating both processes to be important in these disks. However, no 
clear correlation is observed between the mass fractions, suggesting that the two processes 
may proceed independently from each other at different rates (e.g., Bouwman et al. 2001; 
Sargent et al. 2006). In the case of the three Taurus brown dwarfs, CFHT-BD-Tau 8, 
2M04230607 and 2M04400067, that have negligible large grain mass fractions, substantial 
grain growth may have already occurred in these disks, with the larger grains settled towards 
the disk midplane leaving only the small crystalline dust in the optically thin regions (e.g., 
van Boekel et al. 2005). Increasing dust sedimentation would thus result in an increase in 
the fraction of small amorphous grains in the disk surface layers (e.g., Honda et al. 2006; 
DuUemond & Dominik 2008). This is consistent with the >60% small amorphous grain 
fractions found for these objects. These disks may thus be more evolved systems than the 
rest in the sample. The three brown dwarfs however show intermediate flaring, with the 
normalized 13 /im flux densities between 1.03 and 1.18, and thus do not have particularly 
fiat disk structures which would be expected if significant dust settling had occurred in the 
disks. 
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4.3. Dependence on X-ray Luminosity and Disk Masses 

Figures compare the X-ray luminosity witli tlie large-grain and crystalline mass 
fractions. The X-ray luminosities in the 0.3-10 keV range have been obtained from Giidel 
et al. (2007). These observations were part of the XMM-Newton Extended Survey of the 
Taurus Molecular Cloud (XEST) project. Five of these objects have a 'variability flag' of 
0, i.e. they exhibit only low-level variability with no flares visible in their X-ray spectra. 
For the brown dwarf V410 X-ray 6, two measurements are plotted. The measurement from 
XEST exposure number 23, made on March 11, 2001, shows low-level variability (a variability 
flag of 0), and is plotted as a red open diamond. The measurement from XEST exposure 
number 24, made on March 12, 2001, has a flag of 2, i.e. clear flaring is observed in the 
X-ray spectrum (Giidel et al. 2007). The X-ray luminosity is higher by a factor of ~2 for 
this second measurement. The full range in the X-ray emission for these brown dwarfs are 
marked by black open diamonds in Figs. [7^-f. For two other objects, CFHT-CD-Tau 8 and 
2M04141188, the 95% upper limits are plotted. 

Both the large-grain and crystalline mass fractions show some dependence on the strength 
in the X-ray emission. Fig. [7^ shows an increase in the large-grain mass fractions with in- 
creasing X-ray luminosities, while higher crystallinity levels are observed among the weaker 
X-ray sources (Fig. [7f). Several works have concentrated on the effects of X-rays as a source 
of ionisation in protosellar disks (e.g., Glassgold et al. 1997; Ilgner & Nelson 2006). In order 
to explain the angular momentum transport mechanism that drives accretion in protostellar 
disks, Glassgold et al. found that X-ray ionization produces an accreting surface layer that is 
coupled to the magnetic field, and lies over a deeper, quiescent layer. Ilgner & Nelson (2006) 
have considered X-ray irradiation from the central star to be the main source of ionisation, 
and have studied the vertical diffusion of chemical species, that mimics the effects of turbu- 
lent mixing in the disk. These authors have shown that X-ray ionisation induces turbulence 
that results in a decrease or a complete removal of the "dead" zones in the disk. Such zones 
are defined to be the regions that are too neutral and decoupled from the disk magnetic field 
for MHD turbulence to be maintained, as opposed to the "active" zones that are sufficiently 
ionised for the gas to be well-coupled to the magnetic field, and thus able to maintain turbu- 
lence. The main criteria for turbulent mixing to be effective is a significant presence of heavy 
metals (magnesium) in the disk surface layers, since the recombination time between metal 
ions and electrons is longer than the turbulent diffusion timescale. Increasing X-ray irradia- 
tion should thus result in more active turbulent mixing in the disk. This would decrease the 
level of sedimentation and a higher fraction of large-grains would be detected in the upper 
disk layers, as probed by the silicate feature (e.g., Dullemond & Dominik 2008). For the case 
of crystalline silicates, stronger X-ray emission seems to remove crystaline silicates from the 
disk surface layers. V410 X-ray 6, for which even the quiescent emission is stronger than the 
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rest in Fig. [Tf, shows a nearly 0% crystalline mass fraction. However, X-ray flaring may be 
responsible for varying its crystalline structure, as noted recently for the young solar-type 
star EX Lupi (Abraham et al. 2009). A similar anti-correlation between X-ray emission and 
crystalline fractions has been reported by Glauser et al. (2008) among higher mass T Tauri 
stars in Taurus. Amorphization by ion irradiation may be a possible explanation. Such 
processes have been discussed to explain the absence of crystalline sihcates in the ISM (e.g., 
Bringa et al. 2007; Jager et al. 2003; Kemper et al. 2004). Some laboratory experiments 
have shown low-energy (keV) ions to be efficient in amorphizing silicate dust (e.g., Demyk 
et al. 2001; Jager et al. 2003). Though the amophization timescales estimated by some of 
these studies are quite long. In the ISM, Kemper et al. (2004) suggest that the amorphiza- 
tion process occurs on a timescale much shorter than the grain destruction timescale, and 
it would take ~9 Myr to achieve a crystallinity level of 0.4%. Bringa et al. (2007) estimate 
a ~70 Myr timescale for amorphization of ISM silicates by heavy-ion cosmic rays at GeV 
energies. On the other hand, Glauser et al. (2008) suggest that coronal X-ray emission may 
correlate with other high energetic radiation in effectively destroying the crystalline silicates. 
This may reduce the timescale over which such processes take place. 

For three objects in our sample, 2M04141188, CFHT-BD-Tau 6 and 2M04414825, we 
were able to find disk mass measurements of OAdMjup, 0.85Mjup and 0.79Mj;7P, respec- 
tively, from Scholz et al. (2006). The disk masses were derived from the millimeter fiuxes for 
these objects. The estimate for 2M04141188 is a 2-a upper limit. Honda et al. (2006) have 
noted a weak trend (5% significance level) of high crystallinity in objects that are weakly ac- 
creting and have less massive disks. In Fig. [T]!, CFHT-BD-Tau 6 is a weak accretor but has a 
more massive disk, and similar crystalline levels, as the more actively accreting 2M04141188. 
We also do not find a higher abundance of large-grains in the more massive disks; the large- 
grain fractions for CFHT-CD-Tau 6 and 2M04414825 are 23% and 73%, respectively, at a 
similar disk mass of ~0.8Mj[/p. Disk mass measurements for a larger sample are required 
to investigate any such trends. 



4.4. Dependence on Spectral Type 

Figs. [8]a-b show a comparison of the large grain and crystalline mass fractions among the 
Taurus stellar and sub-stellar objects. Data for T Tauri stars in Taurus has been obtained 
from Sargent et al. (2009) and are the "warm" large and crystalline mass fractions, i.e. the 
fraction in the inner regions of the disk at radii around 0.6 AU. We note here that Sargent 
et al. have also considered sub-micron crystalline silicates in their models. Therefore, a 
comparison with their estimated fractions is appropriate. These authors have reported a 
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weak anti-correlation between the large grain mass fraction and the stellar mass, with higher 
abundances for lower mass stars. However, inspecting their Fig. 23 shows a large spread in 
the large-grain mass fractions for stellar masses between ~1 and 0.1 Mq. While we do see 
a similar spread in the large-grain fractions between and ~80% at all spectral types, the 
mean and median fractions for brown dwarfs are found to be 23.4% and 18.2%, respectively, 
and are a factor of ~2 smaller than that found for the earlier type stars (44% and 45%, 
Sargent et al. 2009). The mean and median crystalline mass fractions, on the other hand, 
are found to be comparable among the two sets, 17% and 20% for the brown dwarfs, and 
17% and 11% for the higher mass objects. The spectra shown in Fig. [5^ have not been 
included in calculating the mean and median values. Brown dwarfs are expected to show 
stronger signs of grain growth, as discussed in many previous works (e.g., Kessler-Silacci 
et al. 2006; Apai et al. 2005). As noted above in §4.2^ we have three extreme cases in 
our sample, the objects 2M04230607, 2M04400067 and CFHT-BD-Tau 8, that show less 
than 0.01% large-grain mass fractions, and these may be in a more advanced stage of disk 
evolution than the rest. Excluding these objects results in a mean large-grain mass fraction of 
~30%, which is more comparable to the 44% for higher mass stars. Both the crystallization 
and grain growth processes contribute towards flattening of the observed features in brown 
dwarfs. The flatter features do not necessarily imply a higher degree of grain growth in the 
disks, and thus higher large-grain abundances compared to higher mass stars. In Fig. [HJo, 
no decline in the crystallinity levels is observed towards the later types, as was earlier noted 
by van Boekel et al. (2005) for Herbig Ae/Be stars. The median crystalline mass fractions 
are actually higher by a factor of ~2 in brown dwarfs compared to earlier type stars. As 
discussed in Kessler-Silacci et al. (2007), while the silicate emission zone is larger and lies 
farther away from the central star in higher luminosity objects, the more luminous stars are 
also more capable of heating amorphous silicates located farther away in the disk. Thus 
the crystallization process should proceed independent of the stellar luminosity, resulting in 
similar percentages of crystalline silicates across all spectral types. 



4.5. Age Dependence 

At an age of ~1 Myr, 5 out of the 20 Taurus brown dwarfs studied here show fiat silicate 
spectra, resulting in a 75% detection rate for emission in the 10 /im feature. A comparable 
detection rate of 82% (14/17) was reported recently by Pascucci et al. (2009) for cool star 
disks (SpT ~ M5-M9) in the ~2-3 Myr Cha I star-forming region. Among the older UppSco 
systems (~5 Myr), Scholz et al. (2007) have reported a ~23% (3/13) detection rate. As 
discussed by these authors, brown dwarf disks by this age have undergone significant dust 
evolution, i.e. substantial grain growth followed by dust settling, resulting in featureless 
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silicate spectra for nearly 80% of the disks. Objects such as uscoll2 and some others which 
have flat, very low contrast 10 fim features could then be transitional objects. The above 
comparison suggests a decline in the 10 fim emission detection rate in brown dwarf disks 
with age. However, given the uncertainty on the ages, it is difficult to confirm any such age 
dependence. At an age of ~4 Myr, Sicilia-Aguilar et al. (2007) have found 6/19 low-mass 
stars (SpT G, K, and M) in the Tr37 cluster with marginal emission in the 10 fim feature, 
resulting in a detection rate of ~68%. At a similar age, the detection rate for UppSco brown 
dwarfs is lower by a factor of ~3. This may suggest a shorter timescale over which grain 
growth and dust settling occurs in cool star disks. However, this could again be explained 
by the smaller disk radii probed by the 10 fim feature around brown dwarfs, that would 
show stronger dust processing signatures in these objects compared to higher mass stars 
(Kessler-Silacci et al. 2007; §4.2p . Flat features do not necessarily imply grain growth to 
sizes larger than 10 /im in the disk, followed by sedimentation. Other factors such as strong 
stellar winds could result in a dearth of optically thin material in the upper disk layers, and 
thus fiat silicate spectra (e.g., Bouwman et al. 2001). 

The mean and median crystalline mass fractions for the Taurus brown dwarfs are found 
to be 17% and 20%, respectively. For the Cha I brown dwarfs, the mean and median are 
37% and 46% (Apai et al. 2005). This suggests an increase in crystallinity with age among 
sub-stellar objects. For the UppSco brown dwarf usdl61939, the crystalline mass fraction is 
26%, higher than the Taurus objects. The trend observed among brown dwarfs is contrary 
to the declining trend noted for higher mass T Tauri systems, where the median crystalline 
mass fraction drops from 11% in Taurus to just 2% in the FEPS sample (3-5 Myr; Bouwman 
et al. 2008). The FEPS sample from Bouwman et al. (2008) work constitutes seven pre-main 
sequence systems, and these are the only systems among the 328 targets studied under the 
FEPS Spitzer Legacy program that show signatures of solid-state dust components in their 
mid-infrared spectra. Watson et al. (2009) explain processes such as collisional grinding of 
crystalline grains, or amorphization, that may result in an absence of crystalline silicates at 
older ages. If so, such processes seem to be less prominent in brown dwarf disks. 

The amount of crystalline silicates in a disk is also found to be dependent on the presence 
of dense disk material in the inner disk regions (Sargent et al. 2006). As discussed by these 
authors, if crystalline silicates form due to thermal annealing of amorphous grains in the 
warm inner disk regions, and then radially transported to the outer disk, then the presence 
of an inner hole in the disk would result in a reduction in the crystalline silicate production. 
This may be the case for the brown dwarf V410 X-ray 6 in our sample, that shows a nearly 0% 
crystalline mass fraction. The sub-micron amorphous grain fraction for this object however 
is quite high (85%), and the large grain fraction is only 14%, indicating that significant dust 
processing has not yet occurred in this disk. Among the 65 T Tauri systems in Taurus studied 
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by Sargent et al. (2009) and Watson et al. (2009), 5% are found to have mass fractions 
of amorphous sub-micron grains >95%, with no signs of crystalhne or large sihcates. One 
out of the 20 Taurus brown dwarfs studied here shows a ~87%±10% sub-micron amorphous 
sihcate fraction, resuhing in a similar 5% fraction of the sample still dominated by ISM-like 
grains. At a young age of ~1 Myr, we see a variety in the silicate features in Taurus, from a 
narrow profile with a peak at 9.8 /xm, still dominated by unprocessed dust, to a broad, fiat 
feature peaked at 11.3 fj,m dominated by both large and crystalline grains, to nearly flat, 
low-contrast features. None of the objects in our study show both a neghgible (< 1%) large 
and crystalline mass fraction, and some form of dust processing has occurred in these disks. 
For half of the young brown dwarfs studied here, we find crystalline mass fractions >20%. 
Honda et al. (2006) have noted the crystallization of 5%-20% amorphous silicates in low- 
mass pre-main sequence stars at a very early stage of or before the T Tauri phase. Thus not 
only are such mineral formation processes common among stellar and sub-stellar objects, but 
the timescale over which these take place are also similar, and seem to occur fairly quickly 
and at an early evolutionary stage of the system. And if these processes subsequently lead 
to planet formation, then the building blocks are already present in brown dwarf disks at a 
young age of ~1 Myr. 



5. Summciry 

We have conducted a compositional analysis of the 10 /xm silicate spectra in the Taurus 
and UppSco brown dwarf disks. A comparison of the crystalline and large grain mass frac- 
tions among these two data sets, as well as with brown dwarfs in Cha I and higher mass T 
Tauri stars in Taurus indicates the following: 

1. In comparison with T Tauri and Herbig Ae/Be stars, we find flatter features for the 
brown dwarf disks, consistent with previous findings of heightened fiattening in lower- 
mass objects. 

2. While the strength in the feature shows no correlation with either the large- grain or 
crystalhne mass fractions, the shape is found to be affected by both processes. 

3. The Taurus brown dwarf 2M04141760 shows the least amount of dust processing in its 
disk, with a small amorphous grain mass fraction of ~87%. This is consistent with a 
flared disk geometry observed for this object. 



4. 



For most objects, we find nearly equal fractions of large-grain and crystalline silicates, 
indicating both processes to be important in these disks. 
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5. The large- grain mass fractions are found to increase with an increasing strength in the 
X-ray emission, while the opposite trend is observed for the crystalline mass fractions. 

6. The median large- grain mass fraction for Taurus brown dwarfs is found to be lower 
by a factor of ~2 compared to higher mass T Tauri systems in Taurus. The median 
crystaUine mass fraction on the other hand is found to be a factor of ~2 higher. 

7. No strong correlation is observed between the extent of dust settling in the disk, and 
the large-grain or the crystalline mass fractions. 

8. Among the UppSco brown dwarfs, only usdl61939 has a S/N high enough to properly 
model its silicate spectrum. Our modeling indicates a 74% small amorphous grain 
mass fraction, and a ~26% crystalline mass fraction for this object. 

9. A comparison with the Cha I brown dwarfs suggests an increase in crystallinity with 
age among sub-stellar objects. The detection rate for emission in the 10 /xm silicate 
feature shows a decline with age, with only 20% of the brown dwarfs at ~5 Myr showing 
any detectable emission in the feature. 
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^Reference for spectral type: (1) Luhman 2004, (2) Luhman et al. (2006), (3) Luhman (2006), (4) Ardila et al. (2000), (5) Martin et al. (2004). Uncertainty in spectral type is ±0.25. 

'^Best-fit values for the multiplicative factors in units of lO""*"^ W m~^. 
■^The reduced x^-^alue (dof=651). 

•^Percentage of small amorphous olivine and pyroxene silicates compared to all other silicates. 
^Percentage of large amorphous olivine and pyroxene silicates compared to all other silicates. 
^Percentage of sub-micron crystalline silicates (enstatite and forsterite) compared to all other silicates. 
^Results from remodeling considering large amorphous grains only. 
^Results from remodeling without any crystalline silicates. 
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Table 2. Characteristics of the grain species used for the fitting procedure. 





( '()iiii>ijsi! ion 


Siz(! [//uij 


Sliapi- 


Kclcri-ncc 


/);, [// (7)1 


Amorphous olivine 


MgFcSi04 


0.1 & 2.0 


Homogeneous spheres 


Dorschner ct al. (1995) 


3.71 


Amorphous pyroxene 


Mgo.8Feo.2Si03 


0.1 & 2.0 


Homogeneous spheres 


Dorschner et al. (1995) 


3.71 


Crystalline forsterite 


Mgi.gFeo.iSi04 


0.1 


Irregular (CDE) 


Fabian et al. (2001) 


3.33 


Crystalline cnstatite 


MgSiOs 


0.1 


Irregular (CDE) 


.Jagcr ct al. (1998) 


2.80 


Amorphous silica 


Si02 


0.1 


Irregular (CDE) 


Fabian ct al. (2000) 


2.21 
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Fig. 1. — The normalized continuum-subtracted spectra in units of {F,y — Fc) / Fc. The dashed 
horizontal line represents the continuum. The vertical lines mark the amorphous olivine and 
crystalline forsterite peaks at 9.8 and 11.3 /i.m, respectively. The spectra have been arranged 
in the decreasing order of Fpeak values. 




Fig. 1. — Continued. The second panel shows the UppSco objects. Bottom two panels show 
the spectra for the 5 Taurus objects that do not show any silicate emission. 
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7 8 9 10 11 12 13 14 15 

X(nm) 

Fig. 2. — Normalized spectral profiles for the five dust species considered for modeling. 
Upper panel shows the profiles for small and large amorphous olivine (black) and pyroxene 
(red) silicates. Bottom panel shows the profiles for the crystalline silicates. 
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Fig. 3.— Model-fits for V410 Anon 13 (top) and SST-Lup3-1 (bottom). 
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2M04141760 GM Tau 2M04141188 




Kiim) k{vm) XOrni) 

2M0441 4825 V41 X-ray G CFHT-BD-Tau 9 




Fig. 4. — Model-fit to the lO/xm silicate feature for Taurus brown dwarfs. Colors repre- 
sent the following: red-observed spectrum; black-model fit, grey-small amorphous olivine; 
green-small amorphous pyroxene; cyan-silica; blue-forsterite; pink-enstatite; black dashed- 
large amorphous olivine; orange-large amorphous pyroxene. Thin dashed line represents the 
continuum. The spectra have been arranged in the same order as Fig. [H 




4. — Continued. Bottom panel shows the fits for the UppSco brown dwarfs. 
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Fig. 5. — (a): Model-fits for objects witli negligible large grain mass fractions. Left panel 
shows the best-fit with the lowest reduced-x^ value, obtained using all five dust species. 
Middle panel shows the fits obtained without using any crystalline silicates. Right panel 
shows the fits obtained using only large amorphous olivine and pyroxene grains. 
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Fig. 5. — (b): Model-fits for objects witli negligible large grain mass fractions. Left panel 
shows the best-fit with the lowest reduced-x^ value, obtained using all five dust species. 
Right panel shows the fits obtained using only large amorphous olivine and pyroxene grains. 
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Fig. 6. — Top: (a) The shape vs. strength of the sihcate emission features. Taurus brown 
dwarfs are indicated by filled black squares, UppSco ones by filled red circles. Included for 
comparison are T Tauri stars from Pascucci et al. (2009) and Herbig Ae/Be stars from van 
Boekel et al. (2005). Middle, left (b) Fpeak vs. the large grain mass fraction. Middle, right 
(c) Fpeak vs. the crystalline mass fraction. Blue diamonds indicate Cha I brown dwarfs from 
Apai et al. (2005). Bottom, left (d) Fn.s/Fg.g vs. the large grain mass fraction. Bottom, 
right (e) Fh ^/Fq s vs. the crystalline mass fraction. 
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Fig. 7. — Top, left: (a) Disk geometry for the Taurus and UppSco brown dwarf disks. Also 
plotted are Cha I brown dwarfs (blue diamonds). The values next to the target names are 
the crystalline mass fractions. Top, right (b) Large grain vs. the crystalline mass fractions. 
Middle, left: (c) Disk mass accretion rate vs. the large grain mass fraction. Upper limits are 
indicated by blue arrowheads. The accretion rates have been obtained from MuzeroUe et al. 
(2003; 2005). Middle, right: (d) Disk mass accretion rate vs. the crystalline mass fraction. 
Bottom, left: (e) The X-ray luminosity plotted against the large-grain mass fractions. Red 
open diamond denotes the quiescent state for V410 X-ray 6. Black open diamonds mark the 
range in X-ray emission. Upper limits are indicated by blue arrowheads. Bottom, right: (f) 
X-ray luminosity versus the crystalline mass fractions. 
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Fig. 8. — Left: (a) SpT vs. the large grain mass fractions for T Tauri stars and brown 
dwarfs in Taurus. The value of -2 indicates a SpT of K7, -1 is K5, while 0-9 are M0-M9. 
Right (b) SpT vs. the crystalline mass fractions. 



